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Abstract 

Using the data coming from the new 182 Gold type la supernova samples, the baryon acoustic 
oscihation measurement from the Sloan Digital Sky Survey and the H(z) data, we have performed 
a statistical joint analysis of the DGP brane-world model with a high curvature Gauss-Bonnet 
term in the bulk. Consistent parameters estimations show that the Gauss-Bonnet-Induced Gravity 
model is a viable candidate to explain the observed acceleration of our universe. 
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A variety of cosmological observations suggests a concordant compelling result that our 
universe is undergoing an accelerated expansion, which is one of the deepest theoretical 
problems in cosmology [1]. Within the framework of general relativity, the acceleration 
must be associated with the so called dark energy, whose theoretical nature and origin are 
the source of much debate. Despite the effective negative equation of state uj < —1/3 from 
the robust observational evidence, we know little about the dark energy. 

An alternative approach which does not need dark energy to explain the late-time accel- 
eration is motivated by sting theory via the brane-world scenarios. In the late-time universe, 
one of the simplest extra-dimensional brane-world model which describes the cosmological 
evolution at low energies is the DGP model [2], [3]. In this model, gravity leaks off the 
4-dimensional brane into the 5-dimensional bulk at large scales. Gravity leakage at late- 
times initiates acceleration due to the weakening of gravity on the brane, without the need 
of introducing the mystery dark energy. 

However, the DGP model which modifies Einstein's General Relativity in a consistent 
manner in the infra-red is not free of problems. The most serious one is that such modified 
theories suffer from classical and/or quantum instabilities, at least at the level of linear 
perturbations. Most candidate braneworld models, have been shown to suffer from such 
instabilities or strong coupling or both, [1, 5]. Generically, a ghost mode appears in the 
perturbative spectrum of the theory at the scale where gravity is modified, effectively driving 
the acceleration. Therefore some kind of ultra-violet completion is needed for the DGP model 
in order to be safe at strong coupling. 

There have been some attempts to generalize the DGP model so that they can show ultra- 
violet modifications to General Relativity. One possible way is to introduce a high curvature 
Gauss-Bonnet (GB) term in the gravitational action to display the higher energy stringy 
corrections [6], [7]. An intriguing cosmological model with the combination of infra-red and 
ultra-violet modifications by introducing the GB term in the 5D Minkowski bulk containing 
a Friedmann brane with DGP induced gravity, was presented in [8]. In the general GB 
correction to the Induced Gravity, the late-time self acceleration of the universe is still kept, 
and striking new behaviour in the early universe is also shown [9, 10]. It is of great interest 
to investigate whether such model is a viable cosmological model. 

The pure DGP model was tested using data from various observations [11, 12, 13, 14, 
15, 19]. In this work we are going to impose constraints on the model parameters by using 
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the latest SNIa data compiled by Riess et al [20], the baryon acoustic oscillations (BAO) 
measurement from the large-scale correlation function of the Sloan Digital Sky Survey(SDSS) 
luminous red galaxies [21] in combination with the H(z) data. We will compare our results 
with the cosmological consequences of the DGP model as they were discussed in [13], [14], [15] 
to see the influence of the GB effect on the DGP model and also disclose the value of the 
GB parameter from observations. The GB correction term has been found effective on 
the modiflcation of the cosmological evolution around z ~ 1. This has been reported, for 
example, on the influence of the equation of state either in the modifled RS model or the 
modifled DGP model [16, 17]. 

All the tests we will use to constrain the parameters of our model are for relatively low 
redshift data. It would be interesting to test our model for high redshifts using observational 
data from CMB anisotropics and matter power spectrum. However, this would require the 
knowledge of evolution of density perturbations of our model, a subject which is not fully 
understood even in the pure DGP model [11, 18]. 

Combing the GB term in the bulk with the Induced Gravity on the brane, the Friedmann 
equation on the DGP brane can be described by the dimensionless variables in the form 
of [8] 

A{jh' + l)h^ = (h^ - (1) 

2 

where the dimensionless variables are 7 = h = Hr, fj, = ^p. The conservation equation 
becomes 

/x' + 3/i(l + tu)/i = , (2) 

where the prime denotes d/dr and u is the equation of state. Here r is the crossover length 
scale, which is two times of the value deflned in [13], [14], [15] where = Ml/2M^ and a 
is the GB coupling constant which has the dimensions of [lenghtY'. 

As discussed in [S], the physically relevant self-accelerating solution which is the gener- 
alization of the DGP model exists when < 7 < 1/16 and has the Friedmann equation 



1-87 ^ v/(l-87)^-87^(3/i. + 6) (n ^^^\ 



_ 2167V^ - 367^(1 - 87)(/i. + 2) + (1 - 87)=^ 
' [(l-87)2-2472(/i, + 2)]3/2 • 

In the special case where a = (7 = 0), the (+) branch of the DGP model can be recovered 
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with the Friedmann equation 



\^mH' = -H+^p^. (5) 

7^0+ r 6r 



For the benefit of the following discussion, we rewrite Eq (3) in the form 



, 1-87 ^ v/(l- 87)^-87^(3/1. + 6) / 47r^ 

where /x. = /xo(l + 2;)^. When 2; = we arrive at 



.2 _ 1 - 87 , 87)2 -872(3/10 + 6) / , 47r 



^'^-^^^ 672 co.^.o + yj , (7) 

2167Vg - 367^(1 - 87)(//o + 2) + (1 - 87)^ 
- [(l-87)2-2472(/.o + 2)]3/^ ' 

where /iq = Hqt. Neglecting the GB correction (7 — >■ 0) Eq (6) reduces to 

hm E{zf = + zf + AVLr + ^/mX , (9) 



7^0+ 
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where ilmo = ^fp — p'' A^o = = In our notation which is consistent with that 

used in [S], the crossover factor r = 2rc, which leads to Qr = l/4:Qrc- Then Eq (9) can 
go back to the equation in the pure DGP model described in [13, 14, 15]. Due to the GB 
correction, the cross-over scale obeys [8] 

2Ho^ <r< 4Ho^ , (10) 

while in the DGP(+) limit (7 ^ 0), r ~ 2H^\ It was found that the physically relevant self- 
accelerating solution of the GB correction to the Induced Gravity has a finite temperature 
big bang, since the density /i. is bounded from above [9]. This upper bound is also the 
requirement of real value of the square root in Eq (4) and Eq (6) which is 

, I-I67 + I672 

Requiring 1 0053^^1 < 1, (there is a milder condition: 7 < 1/16 ) we can have 



h - 1 + (12) 

"max — „ ) k 

67 

which is the initial Hubble rate for the model. Then the upper bound density of GB corrected 
Induced Gravity model reads 

/iz < /imax = ^^Lx(2Vl - 127 - 1) < /ima^ , (13) 



where Umax is the initial density with 7 < 1/16. If 7 — ^ 0, DGP model will be restored 
and hmax = f^max = oo. Using the Induced Gravity model with the GB term of the bulk 
to describe the physically relevant self-acceleration, the density upper bound indicates that 
our universe started from a finite redshift z^ax instead of a singularity at z = 00. For the 
universe without dark energy, /i^ = /io(l + z)^ the upper bound on the density leads to 

(l + ^obsf< (1 + 2.^,^3 = ^, (14) 

where Zobs is the redshift we have the observational data and Zmax is the starting moment of 
the universe in this model. 

In the following we are going to constrain this model by using the latest observational 
data, such as the gold SN la data, the BAO measurement from SDSS and combing these 
obsevations with H(z) data. 

The up-to-date gold SN la sample was compiled by Riess et al [ju]. This sample consists 
of 182 data, in which 16 points with 0.46 < z < 1.39 were obtained recently by the Hubble 
Space Telescope (H ST), 47 points with 0.25 < 2; < 0.56 by the first year Supernova Legacy 
Survey (SNLS) and the remaining 119 points are old data. The SN la observation gives the 
distance modulus of SN at the redshift z. The distance modulus is defined as 

^ = m-M = logio Dl + 5 logio + 25 , (15) 

where Dl is the dimensionless luminosity distance and it given hj Dl = {1 + z) 'W^ ' 

From Eq (6) we see that there are two parameters /iq, 7 in the model. Eq (7) tells us 
that /lo is a function of /io and 7. In order to place constraints on the model, we perform 
statistics for the model parameter 

XsNn,f^o,M) = -2 . (16) 

i ' 

The best-fit values of parameter are shown in Table I, where we have done the marginaliza- 
tion of the nuisance parameter M = 51og^o (l^^) + 25. With the best-fit values of /io, 7 
we can get other parameters Hq, a from their relations, which are also listed in Table I. 
In Fig 1(a), we present the contours of 68.3%, 95.4% and 99.7% confidence levels. It is of 
interest to disclose parameters such as Qmo, which have direct physical meanings so that 
we can compare our model with the pure DGP model. Recalling the relation between (/zq, 
7) and {^mo, ^r) and noting that the parameters transformations from (/xq, 7) to {flmo, ^r) 



Table I: 



Test 


7 








ho 


a (in i/g ^unit) 




SNIa 
SNIa+LSS 
SNIa+LSS+H(z) 


U.UUU_Q QQQ 
U.UUU_Q QQQ 


-•^•^^-0.34 

9 OR+0.55 
^■^"-0.42 

^•'-''^-0.33 


Q ir+O.ll 

n 29+°-°^ 

"-'•^^-0.03 

n 97+0.03 
-0.03 


U.UOUZ_Q QQ20 
n nqi 8+0-0029 

n n'^'^'^+0-0024 

U.UOOO_Q QQ25 


9 88+0-10 
^•^^-0.34 

9 oi +0.15 

9 74+0.11 

^•'^-0.10 


n n8fi+0-005 

U.UUU_Q QQQ 
U.UUU_Q QQO 


158.27 
162.92 
174.04 




Figure 1: Probability contours at 68.3%, 95.4% and 99.7% confidence levels for joint parameters: 
(a) (7, /imo)-plane from the gold sample of SNIa data. The point shows the best-fit from the SNIa 
analysis. The red line corresponds to the limit due to Eq (13). (b) For {QmO, r2r)-plane. The upper 
boundary corresponds to 7 = (pure DGP) and the lower boundary is due to Eq (13). 



have non-zero Jacobi determinant "we can obtain the constraint on the physical 

parameters (fimo? ^r) from the SNIa observations. The best-fit values are listed in Table I 
and contours are shown in Fig 1(b). Our analysis shows that if we use only the SNIa data, 
the constrains are not good and the la range is large. 

An efficient way to reduce the degeneracies of the cosmological parameters is to use 
the SNIa data in combination with the BAO measurement from SDSS [21]. The acoustic 
signature in the large scale clustering of galaxies yields additional test of cosmology. Using 
a large sample of 46748 luminous red galaxies covering 386 square degrees out to a redshif 
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Table II: 







Amjn 


SNIa 


'-'•^^-0.03 


159.97 


SNIa+LSS 


'-'•^^-0.02 


162.92 


SNIa+LSS+H(z) 


n 97+0.02 
-0.02 


174.04 



of z = 0.47 from the SDSS, Einstein et al [21] have found the model independent BAO 
measurement which is described by the A parameter 



^BAO 



2/3 



hn 



dz' 

zbao Jo E{z') 

1 /"^'BAO 



^BAO io E{z' 



-| 2/3 

(17) 



The measurement gives A = 0. 469(72^/0. OS)"*^'^^ ±0.017 at ^bao = 0.35. The scalar spectral 
index is taken to be Us = 0.95 through the three-year WMAP data. In our analysis, we 
have investigated the joint statistics with the SN la data and the BAO measuremen. The 
results are shown in Fig 2 (a), (b) where we show the contours of 68.3%, 95.4% and 99.7% 
confidence level for fiQ, 7 and ^Imo, respectively. The fitted parameters with the la errors 
are shown in Table I, where ho and a are obtained from /io and 7. 

It is of interest to include the Hubble parameter data to constrain our model. The 
Hubble parameter depends on the differential age of the universe in terms of the redshift. 
In contrast to standard candle luminosity distance, the Hubble parameter is not integrated 
over. It persists fine structure which is highly degenerated in the luminosity distance [22]. 
Observed values of H(z) [23] can be used to place constraints on the model of the expansion 
history of the universe by minimizing the quantity 

2 / N \Hobsizi) — Hth{zi)f , . 

Xh{1^M=2^ -2 ■ (18) 

The H(z) test on its own cannot provide tight constrain on the model. It is interesting 
to combine the H(z) data with other observational data to obtain tighter constraints on 
the cosmological model. The result on the joint analysis H(z)+SNIa+BAO is shown in 
Fig 3 (a), (b) respectively, la range parameters' spaces are listed in Table I. It is interesting 
to notice that errors of model parameters have been significantly reduced. 




Figure 2: Probability contours at 68.3%, 95.4% and 99.7% confidence levels for joint parameters: 
(a) (7, /imo)-plane from SNIa+BAO data. The point shows the best-fit. (b) For {^ImOi r2r)-plane. 
The upper line is for the limit 7 = while the lower line is due to the limit set by Eq (13). 




Figure 3: Probability contours at 68.3%, 95.4% and 99.7% confidence levels for joint parameters: 
(a) (7, iimo)-p^SLne from SNIa-|-BAO+H(z) data. The point shows the best-fit. (b) For {i^mOj ^r)- 
plane. The upper line is for the limit 7 = while the lower line is due to the limit set by Eq (13). 
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Table III: 



Model 


Test 








ref 


DGP 
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n 1 q+0.06 
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'-'•-'-•^-0.01 




[13] 


DGP 
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"-"•"J -'■-0.06 


n 90+0.03 

'-'■^■-'-0.03 




[14] 


DGP+GB 


SNIa 


^•-■-^-0.04 
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'-'•-'-^-0.01 






DGP+GB 


SNIa+LSS 


'-'•^^-0.03 


n is+0-01 
'-'•-'-•-'-0.01 






DGP+GB 


SNIa+LSS+H(z) 


n 97+0.03 
-0.02 


u-i-J-o.oi 







Table IV: 



Test 


7 


1^0 




SNIa 
SNIa+LSS 
SNIa+LSS+H(z) 


0.0278 
0.005 
0.003 


1 20+°-'^° 

9 9(:;+0.55 
^•^0-0.42 

9 nQ+0.40 
^•^"J-0.33 


1.11 < Zmax < 1-76 
6.86 < Z^ax < 8.06 
10.76 < z^ax < 12.25 



For the sake of comparison, we have also done the fitting to observations by using 
the pure DGP model where the GB correction is absent. Results are shown in Ta- 
ble II. Despite the big difference in the single SNIa data fitting, in the combined analysis 
SNIa+BAO,SNIa+BAO+H(Z), we can see that the GB correction influence the universe 
evolution, however its effect is very small. In Table III, we include the fitting results for 
the pure DGP model obtained in [13], [11], [15]. For comparison, we need to notice that 
Qrc = 4fir- Using our best fit value of Qri^rc), we obtain r = Hq^Qvc^^ = 2.89if(7^, which 
obeys the inequality of (10) due to the GB correction. 

If the value of the parameter 7 is not zero, we can find the maximum redshift at which the 
universe started its existence in our model. Taking the maximum value of 7 from Table I 
and using relations (13) and (14) we calculate the z^ax allowed from observations in our 
model, which we show in Table IV. 

Table IV is actually the result of fitting of our model to the observational data (at late- 
time with small redshifts). The central zero 7 corresponds to the infinite Zmax, where the 
GB effect can be neglected and then the model reduces to the pure DGP. However, what 
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is significant here is that the result is quite sensitive to 7. Even small nonzero 7 will cause 
dramatic change in the evolution of the universe. Table IV takes 7 just at the edge of 
la contour, which illustrates that the effect of GB is significant and quite possible. For 
the modified DGP model by including the GB correction, it was argued that the combined 
induced gravity and GB effects make the universe start at a finite maximum density and 
finite pressure, but with infinite curvature [8]. In other words, the universe described in this 
model does not start at z = 00, but starts at a finite z which is the maximum redshift 
allowed in the model, z^ax in Table IV is this maximum redshift when the combined 
induced gravity and GB effects are taken into account. However, there is a question of the 
high energy-early-time behaviour of the model. The maximum z^ax allowed from Table IV 
is too small to accommodate the conventional CMB formation at high redshifts even if we 
had the technology to calculate such effects in our model (see [21] for such an attempt). 
Therefore, to go to high redshifts region we have to fine-tune 7 to very small values, making 
the contribution of the GB term at late-times negligible. 

In summary, in this work we have preformed a parameter estimation of the Induced Grav- 
ity model with a higher curvature GB term in the bulk proposed in [8]. We have analyzed 
data coming from the most recent SN la sample, LSS observation and H(z) measurement. 
The results show that the DGP model with the GB correction is a viable candidate to ex- 
plain the observed acceleration of our universe. The value of the GB parameter allowed 
by observation is very small giving only small effects to the corrected DGP model. These 
correction effects are sensitive to changes of the GB parameter. A nonzero value of 7 will 
change significantly the cosmological evolution of our universe. However, to make our model 
consistent with the conventional CMB formation at high redshifts the GB parameter has to 
be fine-tuned to very small values. 
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